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u.  A^tTNACT  nAVSTAR  is  a  triservice  satellite  navigation  system  designed  to  pro¬ 
vide  the  means  for  military  aircraft,  ships, and  ground  users  throughout  the  world 
to  obtain  precise  position  and  velocity  information  at  any  time  of  day.  In  any  kind 
of  weather.  Since  the  user's  position  is  measured  by  the  time  delays  between  his 
location  and  several  satellites,  the  additional  time  delay  due  to  tropospheric 
refraction  must  be  taken  into  account.  For  angles  above  5U  the  time  delay,  or 
corresponding  range  error,  can  be  estimated  from  the  local  surface  refractivity, 
if  known. 

The  worldwide  statistics  of  surface  refractivity  analyzed  in  this  report  show- 
that  refractivity  is  correlated  with  site  latitude,  height  above  sea  level,  and 
season,  all  of  which  should  be  readily  known.  Regression  lines  for  range  error 
corrections  based  on  these  parameters  are  derived,  and  the  accuracy  of  the  cor¬ 
rections  are  estimated.  At  5°,  the  minimum  elevation  angle  at  which  NAVSTAR 
is  intended  to  be  used,  the  range  error  i3  about  80  ft/  As  the  elevation  angle  is 
increased,  the  range  error  decreases  approximately  as  i^e  cosecant  of  the  eleva^ 
tion  angle  and  approaches  a  minimum  value  of  about  7  ft  zenith.  Range  error 
corrections  based  on  a  global  average  surface  refractivity  of  324, 8  N  units  have 
a  standard  error  of  about  8  percent  of  the  total  refractive  range  error.  Correc¬ 
tions  based  on  a  regression  model  that  is  a  function  of  site  latitude,  height  above 
sea  level,  season,  and  elevation  angle,  have  a  standard  error  of  6  percent  of  the 
total  range  error.  If  the  surface  refractivity  at  the  user's  location  is  known,  the 
standard  error  of  the  range  error  correction  is  only  about  3, 7  percent  of  the 
total  correction. _ 
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Tropospheric  Range  Error  Corrections 
for  the  NAVSTAR  System 


l.  (introduction 

At  the  request  of  the  Space  and  Missile  Systems  Organization  (SAMSO), 
A?CRL  (LZN)  is  providing  these  worldwide  tropospheric  refractive  range  error 
corrections  for  the  NAVSTAR  system  [formerly  Defense  Navigation  Satellite 
System  (DNSS)]. 

As  projected,  NAVSTAR  will  allow  military  (and  commercial)  aircraft,  ships, 
and  ground  users  throughout  the  world  to  obtain  precise  position  and  velocity  in¬ 
formation  at  any  time  of  day  during  all  kinds  of  weather  conditions.  The  system 
will  consist  of  several  constellations  of  satellites  generating  pseudorandom  noise- 
coded  signals.  The  time  delays  of  the  signals  will  enable  the  user  to  ascertain 
his  range  to  each  satellite,  and  thereby  obtain  sufficient  information  to  determine 
his  position  and  velocity. 

The  accuracy  of  any  system  tends  to  be  limited  by  tropospheric  refraction, 
which  produces  two  main  effects  on  radiowaves:  angular  bending  and  time  delay. 
Angular  bending— deviation  of  the  ray  path  from  a  straight  line— is  due  to  varia¬ 
bility  of  the  refractive  index  with  position  in  the  atmosphere.  The  delay  in  the 
arrival  time  of  a  signal  occurs  primarily  because  the  index  of  refraction  of  the 
atmosphere  is  greater  than  unity,  and  is  to  a  lesser  extent  due  to  lengthening  of 
the  path  by  angular  bending. 

(Received  for  publication  15  April  1974) 
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Thus,  if  the  distance  from  a. ground  position  or  an  aircraft  to  a  satellite  is  to 
be  accurately  determined,  the  time  delay  produced  by  tropospheric  refraction  must 
be  taken  into  account.  This  time  delay  is  often  represented  in  terms  of  an  equi¬ 
valent  range  error  (1  nsec  *  1  ft).  The  tropospheric  range  error  typically  has  a 
maximum  value  of  about  300  ft  for  a  ray  path  whose  initial  elevation  angle  is  0°, 

As  the  elevation  angle  is  increased,  the  range  error  decreases  approximately  with 
the  cosecant  of  the  elevation  angle,  approaching  a  minimum  value  of  about  7  ft  at 
zenith.  At  an  elevation  angle  of  5°,  the  minimum  angle  at  which  NAVSTAR  is  in¬ 
tended  to  be  used,  the  range  error  is  about  80  ft. 

2.  APPROACH 

In  the  most  accurate  of  the  several  methods  used  for  computing  range  error, 
the  error  is  calculated  either  by  ray-tracing  techniques  or  by  numerical  integra¬ 
tion,  based  on  an  actual  refractivity  profile  along  the  ray  path  or  on  a  vertical  pro¬ 
file  at  the  observer's  position.  If  an  actual  profile  is  not  available,  then  a  repre¬ 
sentative  refractivity  model  must  be  resorted  to. 

The  statistics  of  refractivity  as  a  function  of  altitude  for  climatologically 
diverse  regions  show  that  in  the  first  kilometer,  refractivity  decreases  approxi¬ 
mately  linearly;  above  that  altitude,  it  decreases  exponentially.  It  ha3  also  been 
found  that  the  variation  of  refractivity  reaches  a  minimum  at  an  altitude  of  9  km. 
These  findings  form  the  basis  of  the  widely  used  CRPL  Reference  .Atmosphere— 
1958  Refractivity  Model. 1  The  model  consists  of  a  three-section  refractivity  pro¬ 
file;  a  linear  decrease  from  the  surface  to  1  km;  an  exponential  decrease  from  , 

1  km  to  9  km,  at  which  height  it  is  assumed  to  have  a  value  of  105  N  units;  and 
an  exponential  decrease  throughout  the  remainder  of  the  lower  atmosphere.  The 
only  variable  that  must  be  specified  by  the  user  is  the  actual  surface  refractivity 
value  to  which  the  three-piece  profile  is  attached.  Thus,  if  the  surface  refractivity 
for  a  region  is  known,  an  estimate  of  the  range  error  can  be  calculated.  Other 
refractivity  models  have  been  used  but  they  are  for  the  most  part  a  variation  of 
the  CRPL  Reference  Atmosphere. 

3.  RESULTS 

Representative  surface  refractivity  values  were  determined  from  seasonal 
5 -yr  average  surface  refractivity  data  for  268  stations  scattered  throughout  the 

1.  Bean,  B.R.,  and  Dutton,  E.J.  (1966)  Radio  Meteorology,  Dover  Publications, 
New  York. 
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world.  These  were  used  with  CRPL  profiles  to.  compute  range  error  values 'by 
ray-tracing,  Algebraic  expressions  representing  the  variation  of  range  error  as 
a  function  of  surface  refractivity  Ng,  height  above  sea  level  h,  and:  apparent  eleva*? 
tion  angle  j,  -were  derived  for  defining  the  basic  correction  algorithm  for  the 
NAVSTARsystem.  Care  was  taken  to  avoid  using  transcendental  and  trigono¬ 
metric  expressions  that  would  have  added  considerably  to  the  complexity  of  the 
digital  processor  required  in  the  field.  Hence,  the  algebraic  expressions  that  were 
derived  included  only  the  operations  of  addition,  subtraction,  multiplication,  and 
minimal  division. 

3.1  Range  Error  Function 

The  range  error  AR  can  be  expressed  as  a  product  of  three  factor  functions, 
as  follows: 

AR(o,  h,  N  )  =  G(0)  IRh)  Fth,  N  )  (feet)  ,  (1) 

o  9 

where  G(0)  and  H(h)  are  functions  of  only  the  elevation  angle  6  and  the  observer's 
height  h,  and  the  function  F(h,  N_)  is  dependent  only  on  h  and  the  surface  refrac- 
tivity  Ng.  The  factor  functions,  chosen  to  be  of  polynomial  form,  are: 

G<0)  *  <gQ  +  g^"1  +  g20~2  +  g30~3)  [  g4  +  g6(0  -  g5)2]  ,  (2) 

where  0,  in  degrees,  is  valid  only  for  the  range  0  >  5°,  and 

gQ  =  0, 1556  , 

gx  =  130.8926  , 
g2  =  -105.0574  , 
g3  .  31.5070  , 

g4  ».  1.0000  , 

g5  =  30.0000  , 

g6  =  1.0  X10'4; 

H(h>  =  [(bQ  +  b^h  +  B.GSO^+b^h  +  S.eSOl^  +  bgih  +  e.eSO)'3)]  .  (3) 
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where  h  is  -in  thousands  of  feet,  and 


bQ  =  0.00970  . 

b1  ?  -2.08809  , 

b2  =  122.73592  . 

b3  =  -703.82166  ; 


F(h'  Ns)  =  co  I  +  c2  (h+c0>  +  eg 


».-<=«]  [> 


where 


(4) 


c0  =  3.28084  , 

Cj  =  2.07800  , 

c2  =  0.30480  , 

Cg  =  0.00423  , 

c4  =  1,33333  , 

c5  «  1.41723  X 10‘6  , 

c6  =  315.00000  . 

Note  that  although  G,  H,  and  F  contain  a  total  of  seven  implied  divisions,  only 
three  are  actually  required:  (1  /$),  l/(h  +  2.  630);  and  l/(h+c2>.  Once  these  have 
been  carried  out  the  remaining  terms  are  expressed  as  products  of  the  pertinent 
factors,  requiring  only  the  operations  of  multiplication  and  addition, 

3,2  Refractivity  Model 

The  range  errors  caused  by  tropospheric  refraction  can  be  determined  by 
means  of  Eqs.  1  through  4  if  the  NAVSTAR  user  knows  the  elevation  angle  between 
his  position  and  the  satellite,  his  height  above  sea  level,  and  the  local  surface 
refractivity  N  .  If  he  does  not  know  the  local  N  .  he  must  estimate  it. 

A  statistical  analysis  of  our  surface  refractivity  data  showed  that  the  average 
global  surface  refractivity  is  324. 8  N  units  and  that  the  standard  deviation  of 
this  sample  is  30. 1  N  units,  A  further  analysis  of  the  data  showed  that  Ng  cor¬ 
relates  with  site  latitude,  height  above  sea  level,  find  season.  The  following 
regression  line  of  N_  as  a  function  of  these  parameters  can  therefore  be  derived: 

S 

2  2 

Ng(h,  L,  M)  ’  cfQ  -  o^h  -  OgL  +  ofghs  +  ar^Ls  +  or^hc  +  OgLc  , 
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(5) 


where  h  is  In  feet,  L  it  In  degrees,  M  is  the  calendar  month,  and 


«0 

s 

369.03000  , 

*1 

= 

-0.01053  , 

a2 

s 

-0.92442  , 

a3 

= 

0. 00160  , 

°4 

= 

0.19361  , 

a5 

= 

0.00063  , 

a6 

= 

-0.05958  , 

and  E*‘ 

8  =  sin(^M)  , 

c  =  cos  (^jM)  . 

This  expression  describes  the  worldwide  variation  of  Ng  and  has  a  standard  error 
of  ±  17, 0  N  units,  approximately  ±  5  percent  of  the  mean  Ng  value. 

It  was  also  disclosed  that  as  a  factor  in  predicting  N  ,  knowledge  of  just  the 
season  alone  is  as  valuable  as  knowledge  narrowed  to  a  specific  month.  The 
season  was  effectively  used  as  a  variable  in  the  expression  for  the  regression 
line,  with  the  following  assignment  of  M  values. 

Season  M  Value 

W*nter  1. 5 

Spring  4. 5 

Summer  7. 5 

Fall  10. 5 

With  this  model  for  N_,  the  basic  algorithm  for  estimating  tropospheric  range 
s 

errors  ceases  to  be  a  function  of  N_  and  becomes  a  function  of  $,  h,  L,  and  M, 

s 

denoted 

AR (9,  h,  L,  M)  =  G {0)  H(h)  F(h,  L,  M)  ,  (6) 

where 

F(h,  L,  M>  =  c0  +  c2  <h+c0)  +  c3Ns(h,  L,  M)  -  c4]  * 

•  |  1  -  c5  [ns<h,  h,  M)  -g6]  2  |  .  (1) 
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3.3  Accuracy  of  Range  Error  Correction 


The  standard  errors  in  range  error,  calculated  from  the  standard  errors  in 
surface  refractivity  by  using  Eq.  (6),  are  listed  in  Table  1  as  a  function  of  surface 
refractivity  Ng  and  elevation  angle  0.  These  tabulations  represent  only  the  un¬ 
certainties  that  are  due  to  the  statistical  variation  in  surface  refractivity.  Even 
■when  the  surface  refractivity  is  known,  the  model  has  an  inherent  error,  which 
must  also  be  included.  Unfortunately,  the  accuracy  with  which  range  error  can 
be  corrected  on  the  basis  of  the  CRPL  Reference  Atmosphere—lSSS  is  extremely 
difficult  to  determine. 


Table  1.  Uncertainty  in  Range  Error  Due  to  Uncertainty  in  Surface 
Refractivity  (ft/N  units) 


x 

240 

260 

280 

300 

320 

340 

360 

380 

400 

4 

0.197 

0.  188 

0. 178 

0.168 

0,159 

0. 149 

0. 139 

0.130 

0.120 

5 

0.146 

0. 140 

0.135 

0.129 

0.123 

0.117 

0.112 

0.106 

0.100 

6 

0.116 

0.112 

lESEl 

0.104 

0.101 

0.097 

0.093 

0.089 

0.085 

7 

0.096 

0. 093 

:K§l|| 

0.088 

0. 085 

0.082 

0,080 

0. 077 

0.074 

8 

0. 082 

0.080 

0, 078 

0.076 

0.074 

0.072 

0.070 

0.068 

0.066 

9 

0.071 

0.070 

0.068 

0.066 

0,065 

0.064 

0.062 

0.061 

0.059 

10 

0.063 

0.062 

0.061 

0.059 

0,058 

0.057 

0.056 

0.055 

0.054 

12 

0.051 

0.051 

0. 049 

0.048 

0.048 

0.047 

0. 046 

0.045 

14 

0.043 

0. 043 

0.042 

0.042 

0, 041 

0,041 

0.040 

0. 040 

0.039 

16 

0.038 

0. 037 

0. 037 

0.036 

0,036 

0, 036 

0,035 

0. 035 

0.034 

18 

0.033 

0.033 

0.033 

0.032 

0.032 

0.032 

0,031 

0. 031 

0.031 

20 

0. 030 

0.029 

0. 029 

VMM 

0.029 

0.028 

0.028 

0.028 

25 

0.024 

0,024 

0.024 

0.023 

0,023 

0.023 

0. 023 

0.023 

30 

0.020 

0. 020 

0.020 

0.019 

0.019 

0.019 

40 

0.015 

0.015 

0.015 

0.015 

0.015 

0.015 

0.015 

0.015 

50 

0.013 

0.013 

0.013 

0,013 

0.013 

0.013 

0,013 

0,013 

0.013 

GO 

0.011 

0.011 

0.011 

0.011 

0. 011 

0.011 

0.011 

0,011 

0.011 

70 

0. 010 

0.010 

0.010 

0.010 

0.010 

0.010 

0.010 

0.010 

80 

0.010 

0*010 

0.010 

0.010 

0.010 

0.010 

0.010 

0.010 

90 

0.010 

0.010 

0.010 

0,010 

0.010 

0.010 

0.010 

0.010 

0.010 
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Norton  has  calculated  range  errors  for  the  CRPL  Standard  Sample  and 
plotted  the  standard  deviation  of  the  normalized  range  error  AR./R.  for  cases  ■ 
where  the  surface  refractlvity  is  both  known  and  unknown.  If  it  is  assumed  that 
the  standard  deviation  of  the  range  error  for  the  CRPL  Reference  Atmosphere— 
1958  is  the  same  as  for  the  CRPL  Standard  Sample,  then  the  accuracy  of  the  range 

g 

error  correction  for  a  horizontally  stratified  atmosphere  can  be  estimated. 

Errors  due  to  atmospheric  anomalies  su  "h  as  a  nonatratified  atmosphere,  clouds, 
and  precipitation,  must  Also  be  taken  into  account;  these  are  estimated  to  be  a 
maximum  of  about  2  percent  of  the  total  range  error. 

The  standard  errors  of  range  error  corrections  are  listed  in  Table  2  as  a 
function  of  elevation  angle  and  number  of  parameters  available.  It  is  seen  that  if 
Ng  is  known,  the  standard  error  is  approximately  3.7  percent  of  the  total  value 
of  the  range  error  correction.  If  N  is  not  known  but  the  latitude,  height  above 
sea  level,  and  season  are  known,  then  the  standard  error  Increases  to  slightly 
more  than  0  percent.  If  the  latitude  of  the  observer  is  unknown,  the  standard 
error  is  about  7  percent.  If  only  the  global  average  surface  refractlvity  of 
324,8  N  units  is  used,  the  standard  error  increases  to  about  8  percent  of  the 
range  error  correction. 


4.  SUMMARY 

Tropospheric  range  error  corrections  based  on  a  global  average  surface 
refractlvity  of  324.8  N  units  have  a  standard  error  of  about  8  percent  of  the  total 
refractive  range  error.  Since  the  range  error  from  the  surface  to  a  target  above 
the  troposphere  varies  from  about  80  ft  for  an  elevation  angle  of  5°  to  about  7  ft 
for  angles  close  to  zenith,  standard  errors  of  about  6, 4  ft  to  0,  G  ft  would  be 
expected  if  all  systems  distributed  throughout  the  world  used  corrections  based 
on  elevation  angle  only. 

It  has  been  shown  statistically  that  range  error  corrections  based  on  a  regres¬ 
sion  model  that  is  a  function  of  site  latitude,  height  above  sea  level,  season,  and 
elevation  angle,  show  a  standard  error  of  approximately  G  percent  of  the  total 
range  error.  Thus,  a  25  percent  improvement  in  the  accuracy  of  the  ranfje  error 
correction  can  be  obtained  simply  by  using  parameters  that  can  readily  be 
ascertained  for  any  user  location.  The  range  error  correction  can  be  further 


2,  Norton,  K.  A,  (1964)  Effects  of  tropospheric  refraction  in  earth-space  links, 

Proc,  Second  Tropospheric  Refraction  Effects  Technical  Review  Meeting, 
Electronic  SystemB  Division,  TDR-64-103,  1^155-193. 

3.  Altshuler,  E.E,  (1971)  Corrections  for  Tropospheric  Range  Error.  AFCRL- 

71-0419.  - -  - - 


11 


Table  2.  Standard  Deviation  in  Range  Error  Correction*  a*  a  Function 
of  Information  Available  (ft) 


e 

A 

B 

C 

D 

Avg  AR# 

5 

2.96 

5.11 

5.91 

6.75 

83. 8 

0 

2.  60 

4.36 

5.02 

5.72 

'  71.2 

7 

2,30 

3.79 

4.34 

4. 92 

61.8 

8 

2,05 

3,35 

3.83 

4.33 

54.6 

0 

1,80 

3.00 

3.42 

3,86 

48.9 

10 

1.68 

2.70 

3.07 

3.47 

44.2 

12 

1.42 

2.26 

2.  57 

2.90 

37.2 

14 

1.22 

1.94 

2.20 

2.48 

32.2 

16 

1.0? 

1.70 

1.93 

2.18 

28.3 

18 

0. 95 

1.51 

1.72 

1.94 

25.3 

20 

0,86 

1.37 

1.56 

1.75 

22.9 

25 

0. 70 

1.10 

1.25 

1.41 

18.6 

30 

0. 58 

0.93 

1.06 

1.20 

15.8 

40 

0,43 

0.69 

0.79 

0,89 

12.3 

50 

0.37 

0.  GO 

0.68 

0. 77 

10.3 

60 

0.33 

0,51 

0.  59 

0.  67 

9.1 

70 

0.30 

0.48 

0.  54 

0.61 

8.4 

80 

0.29 

0,47 

0.  53 

0.  60 

8.0 

90 

0.29 

0.47 

0.53 

0. 60 

7.9 

'y*, 

0.061 

0.070 

0.080 

A:  N  known 

B:  N„  unknown;  season,  height,  and  latitude  known 

C;  N  unknown;  season  and  height  known 
s 

D:  No  information 
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improved  if  the  surface  refractlvity  at  the  user  location  is  known;  the  standard 
error  of  the  range  error  correction  for  this  case  is  only  about  3. 7  percent  of  the 
total  correction. 

It  should  be  evident  that  very  accurate  range  error  corrections  will  be  required 
only  for  those  locations  where  the  satellites  are  at  low  elevation  angles  since  the 
total  range  error  decreases  as  the  elevation  angle  increases.  Also,  the  range  error 
decreases  as  the  user's  height  above  sea  level  increases  since  the  dense  region 
of  the  troposphere,  which  produces  most  of  the  additional  time  delay,  is  below  the 
user. 

From  Table  2  it  is  apparent  that  the  range  error  is  very  strongly  a  function  of 
elevation  angle.  In  the  work  reported  here  the  assumption  has  been  that  the  user 
knows  the  precise  elevation  angle  between  his  location  and  the  satellite.  For  low 
elevation  angles  this  is  particularly  important  because  the  range  correction  at  those 
angles  changes  rapidly.  How  accurately  the  position  can  be  determined  will  there¬ 
fore  be  very  much  dependent  on  geographic  location. 

It  must  be  emphasized  that  our  results  have  been  directed  toward,  and  are 
mainly  applicable  to,  the  nonstationary  case.  If  the  user  is  confined  to  a  region 
within  which  the  climatology  does  not  vary  appreciably  from  point  to  point,  then 
more  accurate  range  error  corrections  can  be  derived  from  refractivity  statistics 
for  that  region. 
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Errata 


Page  7  -  Change  Equation  (3)  to  read 


H(h)  -  C(bQ  +  b,(h  +  8.6286) +  b2 (h  +  8.6286) "Z  +  b3 (h  +  8.6286) “3)]  , 


Page  8  -  Change  Line  9  to  read 


c.  *  6.81758 


Page  8  -  Change  Line  16  to  read 


three  are  actually  required:  (1/9),  1/ (h+8.6286) ,  and  1/(h+c2).  Once  these  have 
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